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Abstract 

An exact expression for the mass distribution p(M, t) of the ballistic ag- 
gregation model in one dimension is derived in the long time regime. It is 
shown that it obeys scaling p(M,t) = t _4//3 i ? (M/t 2 / 3 ) with a scaling function 
F(z) ~ z~ x l 2 for z 1 and F(z) ~ exp(— z 3 /12) for z 3> 1. Relevance of 
these results to Burgers turbulence is discussed. 

Ballistic aggregation provides a simple model of nonequilibrium statistical physics which 
is a natural version of a dissipative gas of hard spheres where particles follow the basic laws 
of mechanics. It consists in a one-dimensional gas of point-like massive particles which move 
freely until they collide. The perfectly inelastic collision of two masses conserves the total 
mass and momentum, while dissipation occurs as kinetic energy is loss in each collision. One 
can anticipate the formation of more and more massive while slower and slower aggregates. 

This model was introduced by Carnevale, Pomeau and Young [IJ where they conjectured, 
based on scaling arguments and numerical simulations, an asymptotic scaling regime for 
the mass distribution p(M,t) = F(M/(M) t )/(M) 2 t . The average mass per aggregate was 
supposed to grow algebraically with time as (M) t ~ t 2 ^ 3 and the scaling function had a simple 
universal exponential form F(z) = exp(— z) independent of the initial conditions. Later, 
this conjecture was reinforced by Piasecki [[| where he solved the hierarchy of dynamical 



equations governing the system inside a mean-field approximation scheme. 

This system, in its continuous limit, was also studied as a simplified astronomical model 
for the agglomeration of cosmic dust into macroscopic objects [Q. In the ballistic aggregation 
model, the aggregates interact only through their collisions. An aggregation model where 
gravitational interactions are present has been studied in J3J. 

It is important to mention the connection between this model and some solutions of the 
Burgers equation. At very high Reynolds number, the asymptotic solution of the Burgers 
equation consists of a train of shock waves. The laws of motion which govern the dynamics 
of these shock waves are found to be equivalent to a ballistic aggregation system (see ||). 

In this letter, I verify the scaling hypothesis for the mass distribution and find in an 
exact calculation an explicit form for the scaling function. It happens to be different from 
the conjectured simple exponential. 

Rather than solving the set of partial differential equations governing the evolution of 
the system, I exploit the fact that, once the initial state of the system is given, the dynamics 
is completely deterministic. Our approach will thus be based on a statistical study of the 
initial conditions and is largely inspired by the work of Martin and Piasecki |J. 

Initially, particles having all the same mass m are regularly placed on a line with the same 
inter-particle distance a. Initial mass density is thus po — m/a. The initial momentum of the 
thermalized particles are not correlated and are distributed according to the same Gaussian 
distribution 0(p) = yjf3/(2irm) exp(— (3p 2 /(2m)) where I now choose (3 = 1/2 without loss 
of generality. 

I compute now the density distribution p m (X, M, P,t) where p m (X, M, P,t)dM dP dX 
is the number of aggregates located in (X, X + dX) with momentum in (P, P + dP) and 
mass in (M, M + dM) at time t. 

When the coordinates (X, M, P, t) of an aggregate are given, they uniquely define the 
number n = Mj m as well as the initial positions X — Pt/M — Mj (2po) < Xi < X — Pt/M + 
Mj (2p ) (i — 1, . . . , n) of its constituents. A crucial point is that an aggregate, once formed, 
is moving according to the movement of the center of mass (CM) of its constituents, which 
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can be determined from the initial state. I label the location of the CM at time t of the r 
particles located initially at (J + l)a, (j + 2)a, . . . , (j + r)a by 

1 r r + 1 t r 
X r j+l {t) := mx i+i + tPj+i = 3 a + -^~ a + C 1 ) 

The mass distribution can thus be determined from the initial conditions and an aggre- 
gate of mass M = mn is present at time t iff 

• the CM of its leftmost s particles has met the CM of its rightmost n — s particles for 
all s — 1, . . . , n — 1 up to time t, leading to X? +1 (t) > X" t 7/ +1 (t) for 1 > s > n — 1, 

• the CM of the successive groups of particles not constituting the aggregate has not 
met the CM of the aggregate, thus -XJ_ r+1 (t) < X < XJ +n+1 (t) for r > 1. . 

One has (see for details) 



with 9 the Heavyside step function and where M = nm and X = (2j + n + l)a/2 + tP/M. 
The brackets denote the average over the initial distribution of the momentum. 

Using Eq. (0) and the Gaussian form of the initial distribution, one finds the exact scaling 
form 



with M' = M/t 2 / 3 , P' = P/t 1 / 3 and m' = mjt 2 ! 3 . Note that, due to translational invariance, 
the mass distribution does not depend on X. 

Owing to the uncorrelated initial Gaussian distribution of the momentum, one can com- 
pute the density p using an analogy with a Brownian motion in the momentum space under 
particular constraints |||7|], (see Fig. ([!])). One finds 






p m (X,M,P; 



t)=p m (M,P;t) = ^-p m ,(M',P') 



(3) 
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p m ,(M', P') = J m , (-M' - I m ,{M', P')J m , (-M' + (4) 

where J m >(Z) is the probability for a Brownian motion P(r) to start from P(0) = and pass 
above the discrete points P{rm!) > Zrm' — (rm') 2 (r > 1), and I m >(M f , P') is the probability 
for a Brownian motion to start at P(0) = 0, end at P(M') = P' and over-passing the discrete 
points P(rm') > (M' + P'/M')rm' - (rm 1 ) 2 (1 > r > n). 
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FIGURES 
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FIG. 1. The Brownian motion used in the construction of our solution. 

I will derive below an expression for the mass distribution in the limit m' = m/t 2 ^ 3 — > 
which is reached either when t — > oo for a fixed m (asymptotic long time limit) or for any 
fixed time t when m — > (continuous limit). In this limit, one keeps M' and P' of order 
0{1). In terms of the Brownian motion introduced above, the space m! between the discrete 
points barrier shrinks to zero and approaches a continuous barrier which makes the problem 
tractable analytically. Nevertheless, the functions I and J are identically null for m! = 0. 
One should thus keeps track of the first space m! (the Brownian motion will be unrestricted 
up to the first point of the barrier) and will find a mass distribution which is an expansion in 
power of w! . From now on, I drop the subscript ' and set p = 1/2 without loss of generality. 

One finds the dominant contribution in m: 



_ p 2 rco poo 

7 m (M,P)=e-^ / dP l( f ) (P l 

J Mm—m? 



Mm—m 2 



dP 2 <t>{P 2 )K M {P u m, P 2 , M-m) 



m 



d 2 



= P 2M 

7i dP l dP 2 



K M (m,P u M-m,P 2 ) 



+ 0(m 2 



(5) 



Pi=P 2 =0 



and 



_ roc poo 

J m {Y)= dP^Pt) lim/ dP 2 K Y (P h m,P 2 ,Nm) 

JYm-m 2 N^oo JYNm-(Nrn) 2 



'Ym-m 2 N-*oo JYNm-{Nm) 2 

-lim / dP 2 —K Y (P l7 m,P 2l Nm) 

7T W— >oo JYNm-(Nm) 2 6P\ 



+ C(m 3 / 2 ) (6) 

P 1= 



where Kz(Pi, T\, P 2 , t 2 ) is the probability for a Brownian motion to start at P{t\) = P\, end 
at P(t 2 ) = P 2 while staying above the continuous barrier P{t) > /(r) = Zt — r 2 . 
Defining the stochastic process Qij) = P{r) — f(r), one has || 

K z (Pi, n, P 2 , r 2 ) = G(Q U n , Q 2 , r 2 ) exp Q(Qi/'(ti) - Q 2 f(r 2 ) - I jT rfrf(r) 2 )) (7) 

where Qi = Pi — /(t$), (i = 1, 2). The function G is the solution of the equation |J 

- - f GW^n.Q,,^) = (8, 

with G(Qi, r, Q 2 , r) = 5{Qi - Q 2 ) and G(0, n, Q2, ^2) = G(<5i, n, 0, r 2 ) = 0. In our problem 
Ht) = 2. 

The equation (||) can be solved (see [0] for details) and one finds 

G(Q 1 ,t 1 ,Q 2 ,t 2 ) = 2^e fcl2 ^ -— — — (9) 

k>l l A1 \ U k)) 

where Ai is the Airy function || which has an infinite countable numbers of zeroes — on 
the negative real axe (— iO\ = —2.33811 . . . , — uj 2 — 4.08795 ...,...). This function had to 
be expected in this problem as it is known that it arises in the description of a Brownian 
motion with a parabolic drift ||. 
Using Eqs.(|||7|,|9]), one gets 

I m (M) = -e- A/3 / 12 J(M) + 0(m 2 ) (10) 
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with 



Z(M) = 5>-^ M (11) 
fe>i 



In the same way, I use Eqs.(g,0,g) and obtain 

JY/2~M) a /3-(Y/2)' i /3 f°° j„ -x(Y/2-M) „-UiM ~ U k) 



j m{Y ) = J- lim e W2^) 3 /3-(W/ 3 f 00 d xi Y/2-M) e~^ M A X Wk J + 0(m 3 / 2 ). 

(12) 



Using the integral representation of the sum in (0) ||, 

^ M M(?-uJ = r+™ M Ai(z + x) 

jfe Ai'(-cu fc ) Jc-ioo Ai(z) 1 ; 

with c > — u\, one can exchange the integration and find after a tedious analytical calculation 



J m (Y) = ™ e -™ 3 / 3 J(Y) + G(m 3 / 2 ) (14) 

V 7T 



with 

where c > — 

Now, inserting the expression for I and J in Eq.([|), one has in the original variables 
m 2 JM\ / M P^/ 3 \ / M Pt x l z \ (m 2 \ 

p - *) = wv 1 {w) J [~w + ) J \-w ~ j + l [w j («o 

From this equation one has that the concentration c(t) of aggregates, the aggregates 
average mass and momentum and the mean energy per unit of length E(t) behave, for time 
t 1, as 



c(t)~t~ 2 / 3 , (M) t ~t 2 '\ (P) t = 0, J(P^) t ~t l l\ E{t)~r 2 /\ (17) 

A careful integration over P [0] leads to the mass distribution, which is the main result 
of this letter, 

*')4 f © +, B < 18) 

where one sees that it obeys the expected scaling form with a scaling function 
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777 

F(M') = 2—M'l(M')TC(M') (19) 

IT 2 



where 



I i-c+ioo e~ M ' z 



k>l 



X (M <) = £e-"', H{ , n = _J^ d (20) 



with c > — u\. 

The scaling function F(M') is plotted on Fig.0. 



0.20 rr 




FIG. 2. The rescaled mass distribution t 4 / 3 p(M, t) function of M/t 2/3 . 

One can compare the obtained scaling function with the conjectured one (F con j (M') = 
exp(— M')) 0]. In particular small and large arguments present strong differences. Indeed, 
for M' <C 1, one get Ti(M') = 1 + 0(M') while one can estimate X(M') using the asymptotic 
properties of the zeroes of the Airy function Uk = {(3nk)/2\ 2 / 3 + 0(k~ l l 3 ) and find X(M') ~ 
(2 V ^FM' 3 / 2 )- 1 . One thus find 

777^ 1 

f(m') = - kj - 2 - 7 = + o(Vm^), (M'«l). (21) 

One can conclude, for example, that the number N(M ,t) = J M ° p(M, t) dM of aggregates 
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of small masses M < M <C t 2 / 3 at time t is well underestimated by the conjectured form 
which leads to N(M ,t) ~ M /t, while the exact solution gives N(M ,t) ~ y/Mo/t. 

For M' ^> 1, one can estimate the function 7i(M') by the steepest descent method and 
find H(M') ~ v/7rAf /3 / 2 exp(-M /3 /12). On the other end, one has J(M') ~ exp(-wiz) and 
finally 

F(M') = -l-M' 5 / 2 e- M '^- M ' 3 / 12 , (M' > 1). (22) 

7T 3 ' 2 

This is again different from the conjectured function as large masses M ^> t 2 ! 3 have a much 
smaller chance to be present in the system. 

Notice that the asymptotic behaviors of the scaling function are compatible with the 



exact bounds found for the burgers problem |]T0[ with white noise initial condition. On 
the other end, Eq.flTSD solves the shock strength distribution questioned in [B| and studied 
numerically in |TT] . 



It is instructive to compute, along the same line, the collision frequency h> 2 (Mi, M 2 ,t) 
where i>2 dM\ dM 2 dt is the number of collision per unit of volume between masses in 
(Mi, Mi + dMx) and (M 2 , M 2 + dM 2 ) in a time in (t, t + dt). I find 



/OO POO 
dPt / dP 2 
-oo J ~ oo 



Mi M 2 



^(M7-t>- M - M2 



xt~ 2/3 J m > (--^7 - M[j l m \M[,P[)I ml {M' 2 ,P' 2 )J m , U± - M^j , (23) 
with Ml = Mi/t 2 / 3 , P( = P l /t 1 / 3 and vnl = mjt 2 ! 3 , leading to 

(24) 

with 

T{M[, M' 2 ) = (M[ + M^)M(M^J(M{)X(M2)H(M( + M 2 ) (25) 

and X and H, as above. This collision frequency clearly does not factorize in a product of 
functions of Mi and M 2 , respectively. This fact invalidates the assumption on which the 
mass distribution was computed in 0. 
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One can inquire about the universality of these results with respect to other initial 
conditions. Let us first consider a Poissonian distribution of the particles initial positions 
with an average interparticle distance a. The discrete points over which the Brownian motion 
should pass in the construction of our solution are still distributed on the same parabola 
but with irregular spacing. In the long time limit and after rescaling, the spacing between 
points of average a' = a/t 2 ^ 3 become smaller and smaller up to be, in first order in m', 
a continuum. The difference between irregular and regular spacing is thus asymptotically 
erased and the result Eq. ([l8|) should be recovered in this case. 

A bimodal momentum distribution (p{p) = (5(p — po) + 5{p + po))/2 is used in the initial 
state in jjj]. I believe that this should not affect the form of the mass distribution ([T8|) as 
the random walk initiated by this distribution is well approximated, in the long time limit, 
by the considered Brownian motion. 

One can define a distribution where momentum are initially correlated. In this case, one 



expects the scaling function to be different, at least for small M' [ [T^j . 

In summary, I found an exact asymptotic solution for the mass distribution of the ballistic 
aggregation in one dimension. Such an exact solution is not frequent in a nonequilibrium 
system and has permitted to verify scaling hypothesis for this system. While the average 
mass per aggregate was proved to behave with time as (M) t ~ t 2 ! 3 for t 3> 1, as expected 
from previous studies, the scaling function is shown here to be different from the conjectured 
one. This distribution also solve the shock strength distribution of the one dimensional 
Burgers equation in the inviscid limit with a white noise initial condition. 

I gratefully acknowledge numerous useful discussions with P. Martin and J. Piasecki and 
financial support from the Swiss National Foundation. 
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